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Topological aspects of superconductivity in quantum spin-Hall systems (QSHSs) such as thin 
layers of three-dimensional topological insulators (3D Tis) or two-dimensional Tis are in the focus 
of current research. We examine hybrid QSHS/superconductor structures in an external magnetic 
field and predict a gapless superconducting state with protected edge modes. It originates entirely 
from the orbital magnetic-field effect caused by the locking of the electron spin to the momentum 
of the superconducting condensate flow. We show that such spin-momentum locking can generate a 
giant orbital g-factor of order of several hundreds, allowing one to achieve significant spin polarization 
in the QSHS in the fields well below the critical field of the superconducting material. We propose a 
three-terminal setup in which the spin-polarized edge superconductivity can be probed by Andreev 
reflection, leading to unusual transport characteristics: a non-monotonic excess current and a zero- 
bias conductance splitting in the absence of the Zeeman interaction. 

PACS numbers: 72.25.Dc, 73.23.Ad, 74.45.-l-c 


Introduction. Spin-Hall effects are one of the most ac¬ 
tive fields in modern solid state physics m- In partic¬ 
ular, the quantum spin-Hall effect [3 El 0 Ho] allows one 
to generate and convert charge and spin currents in pro¬ 
tected edge channels m- Combining quantum spin-Hall 
systems (QSHSs) with superconductors (SCs) leads to a 
broader spectrum of interesting observable phenomena 
[I2HT4|. These include quantum interference effects re¬ 
ported in mini, indicating superconducting transport 
through the edge states in the QSH regime. Understand¬ 
ing edge superconductivity in QSHS/SC hybrids is also 
instrumental to the proposals to realize Majorana zero 
modes in topological insulators (see, e.g., m and re¬ 
views [niHT])- 

In this paper we predict a unique magnetic-field re¬ 
sponse of QSHS/SC hybrids which is characterized by 
very large effective g-factors reaching the order of several 
hundreds. It originates from the locking of the electron 
spin to the momentum of the superconducting conden¬ 
sate flow generated by an external magnetic field. We 
show that this orbital effect has the form similar to the 
Zeeman spin splitting in thin superconducting films |18j , 
but involves an effective g-factor determined by the pa¬ 
rameters of the QSHS/SC structure, viz.: the edge-state 
velocity, v, the thickness of the SC material, dsCi ^nd the 
London penetration depth, A^: 


9 * = 


2mv 

—-—Xp tanh 
n 


dsc 

2Xl’ 


( 1 ) 


where m is the electron rest mass. Equation 0 sug¬ 
gests a simple way to engineer through an appro¬ 
priate material and structure choice. For typical lat¬ 
eral HgTe- or Bi 2 (Te,Se) 3 -based/Nb systems we predict 
g* « 200-7400. Such giant effective g-factors allow one to 
control the proximity-induced edge superconductivity by 
weak magnetic fields that do not destroy the SC order 
parameter or violate the protection of the edge states 


against backscattering. In the following, we illustrate 
these points in detail by developing a microscopic theory 
for the orbital spin splitting of the superconducting edge 
states and spin-dependent Andreev reflection. 

Model. We consider a lateral QSHS/SC structure (see 
Fig. 0 with proximity-induced superconductivity in the 
QSHS described by the Bogoliubov-de Gennes (BdC) 
Hamiltonian 


^BdG 


/ LIn iAroCTyN 
\-iATocry -H;^ ) ’ 


( 2 ) 


where A is the induced pairing potential, and is 
the Hamiltonian of the normal QSHS. For the sake of 
concreteness, we assume that the QSHS is realized in a 
thin film of a 3D Ti [T9| and described by the four-band 
Hamiltonian ^D] : 


= v{px + eA3,)T^(T3, -I- vpyT^ay - /TToCTo + TtxCJo- (3) 

It models top and bottom surface Dirac states hybridized 
in a thin slab with the coupling energy T. We use the 
representation in which the hybridization term is propor¬ 
tional to the off-diagonal Pauli matrix acting in the 



FIG. 1: Schematic of a QSHS/SC structure on a substrate. 
The system occupies semispace —oo < x < oo,t) < y < oo. 
Inset: QSHS realized in a thin-film topological insulator with 
hybridized top and bottom surfaces (see also text). 
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basis of top and bottom states, while the Dirac terms 
are diagonal in this basis and proportional to the Pauli 
matrix r^. Pauli matrices, ax and ay, represent the spin 
of the surface carrier (cto is the corresponding unit ma¬ 
trix), p = —iH[dx,dy,0] is its momentum operator, and 
p, is the chemical potential. Additionally, we introduce 
an external magnetic field B = [0, i?y,0] and the vec¬ 
tor potential A = [^ 3 ,( 2 :), 0, 0] related by By = dzAx{z) 
and choose e > 0. The function Ax{z) must be deter¬ 
mined by taking into account the screening of the ex¬ 
ternal field by the diamagnetic (London) current in the 
SC [H]. Since realistically the thickness of the QSHS, 
dN, is much smaller than dsci the vector potential in 
the QSHS can be approximated by the Ax value taken 
at the SC surface Ax « —tanh(dsc/2AL) [U]. In 
this approximation, the effect of the magnetic field is 
equivalent to a linearly varying superconducting phase, 
4>{x) = {2e/h)Ax X, whose gradient generates a conden¬ 
sate flow with the momentum 2eAx. 

The relevant symmetry operations for our system are 
time reversal (TR) f, particle-hole (PH) conjugation C, 
and their combination tC. Introducing Pauli matrices 
in particle-hole space T^x,y,z and the complex conjugation 
operation K, we can define the symmetry operations and 


their action as 



f = -iPoT^ayk, HBdoiA) 

HBdG{ — A), 

(4) 

c = -TTxToaok, HBdoiA) 

—HBdciA), 

(5) 

HBdoiA) 

—EBdci—A). 

(6) 


For A ^ 0 the TR is broken, while the PH symmetry is 
preserved, so that the excitation spectrum remains sym¬ 
metric around the Fermi level. 

Orbital spin splitting. We seek the solution to the BdG 
equation iFedC i’ix, y) = E-iplx, y) in the form of a plane 
wave li^ix, y) cx propagating along x with a wave 

vector k and decaying exponentially for ?/ > 0 on a scale 
6. From the requirement that the normal component of 
the current vanishes at the edge, jy{x,y = 0 ) = 0 , we 
derive the boundary condition 

tp{x,y = 0) =TToTyayilj{x,y = 0). (7) 

It is specific to Dirac fermions and preserves TR. In 
we prove that the edge solution exists only for an inverted 
band gap F < 0, with 6 = hv/\r\, and has the following 
structure 

'(p{x,y) = {XsAcryXs,^sAcry^sf ( 8 ) 

We took advantage of the one-dimensional character of 
the edge spin-momentum locking and used the eigen¬ 
states Xs = (Sj 1 )^ of the spin matrix ax with the eigen¬ 
values s = ±1 as well as the related hole spinors (see 




(c) 


(d) 


FIG. 2: Edge dispersion for different parameters. Solid and 
dashed curves correspond to the opposite-spin states with 
s = ±1, respectively. Red and blue colors schematically in¬ 
dicate particle and hole branches corresponding to those in 
the normal system. For By = 0 there is an excitation gap 
at two Fermi points ±fc_F indicated by the vertical lines. For 
By ^ 0 the gap is reduced at each Fermi point (b), vanishing 
completely for By = B, (c). For By > B, (d) the spectrum is 
gapless, albeit each spin branch separately has a gap at finite 
energies. 


EH). The edge dispersion is 

{k) = sveAx ± \/+ (p — shvkY 

= —sg^,pBBy ± A'^ + {p — shvky. (9) 

We note that the spin-momentum locking couples the 
spin s to the condensate momentum cx Ax, giving rise 
to the orbital spin splitting of the levels with the ef¬ 
fective g-factor g^ introduced earlier in Eq. 0 [yB 
is the Bohr magneton]. Let us estimate g* for Bi 2 X 3 
(where X=Se, Te) and HgTe-based structures which typ¬ 
ically have hv « 300 meV nm [HI |H] and, consequently, 
g* « 7.9 nm“^ x Al tanh(dsc/2AL). For commonly used 
SC Nb with Al ~ 50 nm [21] and nanoscale thicknesses 
dsc = 50 -b 250 nm, we obtain 5 * « 200 -b 400, which 
exceeds the electron spin g-factor by at least two orders 
of magnitude. In contrast to the Zeeman spin splitting 
[IS] , the orbital effect considered here does not require 
high magnetic helds to generate spin splitting, as will 
be shown below. At ^ = A = 0 Eq. Q de¬ 

scribes a spin-degenerated Dirac cone. A finite chemi¬ 
cal potential shifts the electron and hole branches rela¬ 
tive to each other, so that the SC pairing opens a gap 


at two Fermi points Akp = Ap/{hv), cf. Fig. 2a for 
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Ax = By = 0. Note that the Fermi points -iikp corre¬ 
spond to the opposite-spin states with s = ±1. There¬ 
fore, the magnetic field splits the states at different Fermi 
points, shifting them relative to each other in energy, 
as depicted in Fig. When the orbital energy scale 
matches the induced gap at \Ax\ = A/(ew) or, equiva¬ 
lently \By \ = S*, with 


i?* = 


cvXl tanh(dsc/2AL) 


eCscAi tanh(dsc/2AL) ’ 

( 10 ) 


where ^sc is the coherence length, the excitation gap van¬ 
ishes [Fig. 2c . For the induced gap A ^ 0.1 meV pK] 
and Nb-structure parameters used above, ^ 5.5 mT 
which is much smaller than the critical field of Nb. For 
By > B^ the spectrum remains gapless. Note, however, 
that each spin branch separately has a gap at finite en¬ 
ergies [Fig. 1^, which is the origin of spin-dependent 
Andreev reflection discussed later. 

Symmetries and protection. Let us examine the prop¬ 
erties of the edge solutions at fixed energy E. For that 
purpose, we solve ([^ for the wave vector 

= skp + j3\/ {E — sveAx)^ — A‘^/Hv, 13 = ±1, (11) 
and find the energy-dependent eigenfunctions as 


a , . ^ exppfc^^x - y/S] f 
^EsVy) fET \ /5 ^ 

^/25 


( 12 ) 


It is expressed in terms of the normal-state edge solution 
<f)s = (s, 1,1, —s)"^ and the BCS coherence factors 


A 




'^Es = s sign[£)ss]u^f, Des = E - sveAx- (13) 

Using k^^{Ax) = -k^^^i-Ax), u§^(A^) = u~/_^[-Ax) 
as well as v^g(Ax) = —v'^^^{—Ax) the action of TR on 
the eigenstates is described by 


'0Bs(^a;) 




(14) 


For Ax = 0 the states V’ss ^'^Es Kramers’ part 


ners. Using (131 the corresponding relations for PH and 


tC are derived: 


ijj' 


P c 

Es ^ 


-s signDEstp_E-s 


i’Es ^ signDEs exp 


. 2/3Vj4L-A^ ' 
hv 


(15) 

TTzToCToV'Ef ■ 

(16) 


In contrast to (14) and (15) the relation for tC is not 


broken by a finite vector potential and connects states 


at the same energy. The price is an x-dependent phase 
factor and a transformation in particle-hole space. C and 
tC both change the quasi-particle character. 

In the following we discuss the protection of the supercon¬ 
ducting edge states against elastic scattering, assuming a 
non-magnetic disorder potential iLdis = U(a:) 7 r 2 rocro- To 
that end, we calculate the matrix element of iLdis between 
states with the same energy E: 


Ves- 


Bfiis 




Es2 


V{x) 


■ exp 


KkEs2 - kEs,)x 


Pi 


•(1 + S 1 S 2 ) ^^ 2 ) ■ 


(17) 


The matrix element vanishes for the states with opposite 
spin projections si = —S 2 - This means the absence of 
backscattering that couples different Fermi points kp = 
EpL/hv. This is true irrespective of the presence or ab¬ 
sence of the magnetic field. Furthermore, there is no 
scattering between the states near each Fermi point, kp 
or kp. For such states Si = S 2 = s, /3i = —132 = /3, and 
{ues)*''^Es ~ ^'^EsY'^'e^s ^ yielding vanishing matrix 
element This can be interpreted as protection by 

the tC symmetry, which has the same matrix structure 
as iLdis, leading to a generalized Kramers’ theorem 






(18) 


Physically, the fC protection means that the states with 
opposite /3s (i.e. particles and holes) cannot be converted 
into each other without being Andreev reflected [26] . For 
By > B^ ( |l 0 | there are four protected zero energy states. 
At the transition By = B^ their number is reduced to 
two: 


, 0 . , exp[iskFX-y/S] / tj), 

{-signAxcl^.s 


(19) 


These are equivalent to two Majorana zero modes. In¬ 
deed, defining 77 = (I -|- signAa;)/2, one can construct 
new states 71 = exp[i 777 r/ 2 ]('(/;° + V'-s) 72 = 

—iexp[i 777 r/ 2 ](' 0 ° — V:®,,) with self-adjoint properties of 
emergent Majorana fermions. 

Proposed experimental detection. The spin split¬ 
ting in the QSHS can be characterized by the spin- 
dependent density of states (sDOS) defined as Ps{E) = 
-/o°°d?/iTr [lmGf(?/,y)]^^ where [imGf (y, y)] is 

the imaginary part of the quasiparticle retarded Green’s 
function |^. Pe{E) diverges at the band edges E = 
sevAx T A and approaches \/{hv) away from the gap. 
The sum p^(U) -|- p\^{E) yields the DOS, while the dif¬ 
ference p^{E) — p^{E) characterizes the spectral spin po¬ 
larization. These quantities are plotted as a function of 
E and By in Fig. The gray shaded planes indicate 
By/B^ = 0, 1 and 1.5 used in Fig. At zero mag¬ 
netic field the DOS shows a quasi-particle gap between 
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FIG. 3: The DOS (a) and spectral spin polarization (b) as a 
function of By and E. The gap closing is clearly visible at 
By = Bf The shaded planes indicate the By values used in 
Figs |2a| and |2d| respectively. 


E = ±A. Within the bands the DOS is constant and spin 
degenerated. Hence the spin polarization vanishes. A fi¬ 
nite By shifts the energy of the s-branches by —sg^UBBy. 
This is reflected by a splitting of the peaks in the DOS. 
At gap closing (101 the two peaks cross and the gaps for 
the different s-branches separate in energy. In the en¬ 
ergy range, where only one s-branch is gaped, a finite 
spin polarization arises. 

Although the spin s, which characterizes the states, 
is derived from a ax eigenstate Xsi the structure of the 
solution Q always combines parts with opposite ax ex¬ 
pectation values. The reason is that the eigenstates of 
^ consist of a perfect mixture of the two surface states 
with opposite helicity. Hence it is to be understood as a 
pseudo-spin and the spin polarization does not result in 
a spin signal, which is directly measurable by a magnetic 
tip, like it was found for HgTe quantum wells HU. How¬ 
ever, comparing Figs.|3a1and|3b|we see that the sDOS is 
identically mapped on the regions with unit DOS, being 
therefore an indirect measure of the spin polarization. 
A direct measurement of the spin polarization can be 
realized in transport experiments. We propose the Y- 
forked 3 terminal normal-superconducting (NSC) junc¬ 
tion, sketched in Fig. [i^ If the system is large enough, 
the counter-propagating edge states with opposite spin 
do not hybridize and can be treated independently, while 
perfect Andreev reflection occurs at the NSC interface 
EauH]. When a bias V is applied across the NSC junc¬ 
tion we calculate the excess current Jex as well as the 
dl/dV characteristics within the BTK formalism [2T1 [26] 
at finite temperature k^T — 0.1 A. Since the spin orien¬ 
tations are in one-to-one correspondence to the geometri¬ 
cal edges, the two spin branches can be detected directly 
by measuring contacts 1 and 2 independently. In Fig. 
we plot /ex as a function of By. While the total cur¬ 
rent (blue dotted line) is symmetric, the individual spin 
branches are odd and show a distinct, spin-dependent 
maximum at By = —seV/{2g^fiB) |21) . Further we show 



/(=2A, eV=2A,/tj,T=0.1A 



By [B.] 


(a) 


(b) 


/i=2A, eV=0.5A, *sT=0.1A ^=2A, eV=1.5A, ArBT=0.1A 



FIG. 4: (a) Y-shaped junction between a normal (N) and 
a superconducting (SG) 3D Ti thin film. The different spin 
branches, depicted as black, dashed (s=l) and red, solid (s=- 
1) lines in (b),(c) and (d), can be measured independently in 
the leads 1 or 2 and summed up to the total signal (blue, dot¬ 
ted). The excess current 7ex (b) is odd in By for a single spin, 
showing a clear maximum at finite By. The d//dF charac¬ 
teristics have a pronounced maximum at zero field, which is 
split when eV > A (d). 


the d //dV characteristics as a function of By in Figs 4c 


and |4d| The spin-dependent quasi-particle gaps corre¬ 
spond to the conductance plateaus and can be detected 
directly by measuring the non-local conductances in lead 
1 and 2 separately. The peak in the total conductance 
at zero field splits when eV exceeds the superconducting 
gap, i.e. when the quasi-particle gap is closed. This is a 
rather general signature of the spin-split superconducting 
states (29) . Since energy and magnetic field enter the so¬ 
lutions (131 in the same way, an analogous signal can be 
obtained by varying V for fixed By. The non-monotonic 
excess current and the split conductance peak are both 
hallmarks of the spin polarization in a superconducting 
QSHS. 


In conclusion we have demonstrated the existence of 
a hybrid superconducting quantum spin-Hall system, 
which is robust against elastic backscattering in finite 
magnetic fields. With appropriately chosen structure 
and material parameters, the system is characterized by 
very large effective g-factors reaching the order of several 
hundreds, allowing one to achieve spin polarization by 
applying a rather weak external magnetic field. The he¬ 
licity and the spin polarization of the emerging states can 
be experimentally detected in 3-terminal NSC-junctions. 
We predict a non-monotonic behavior of the excess cur- 
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rent as well as splitting of the zero field and the bias 
differential conductance peaks. Both features are con¬ 
nected with the closure of the quasi-particle gap at very 
small magnetic fields. 
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for the financial support the German Science Founda¬ 
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FIG. 5: Geometry of SC/QSHS hybrid structure. The dashed curve shows the position dependence of the vector potential in 
Eq. ([2^. 
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Details of the model 


The system is described by the Hamiltonian i?BdG given by Eqs. ([^ and ^ . We introduce three independent sets 
of Pauli matrices, tt^, Tj, and ai, to represent the particle-hole, top-bottom, and spin degrees of freedom, respectively, 
and write H^dG as 

HsdG = hvkxTToT^ax + veAxTTzTzCrx + hvkyiizTzay + veAyiToTzay - p^nzToao + TnzTxao - AiTyToay. (20) 

In this subsection we determine the vector potential A^,, taking into account the screening of the external field by the 
diamagnetic current in the superconductor (SC). We will also show that the vector potential acting on the states in 
the quantum spin-Hall system (QSHS) can be approximated by the value of A^ at the SC surface. 

We choose the z = 0 plane at the SC/QSHS interface (see also Fig. and assume that an external magnetic field 
is applied along the y axis. It is convenient to use the London gauge 

A = Az;{z)x, (21) 


in which the SC order parameter can be chosen real, so that the screening current density can be described by the 
London equation 


3 = 



( 22 ) 


where m* is the effective electron mass, and Ug is related to the density of the Cooper pairs in the SC. Then, Ampere’s 
law in our geometry yields 


n 1 

dlAz,{z) = ^o^A,(^) = y^A,(z), 
m 


(23) 


with the London length \l = \/rn*/[ybdUge^). The boundary conditions to Eq. (23) are symmetric with respect to 
the applied field; 


dzAz,{z = t)) = By, dzAz:{z = dsc) = By. 


(24) 


Here dsc is the thickness of the superconductor in ^-direction. Consequently, the solution of this boundary problem 
is anti-symmetric with respect to the middle of the SC: 


Ag;{z) = By\L- 


sinh 

w( 

z-<^) 

cosh 

dsc 

2\l 



0 < z < dsc- 


(25) 

















7 


By continuity, the vector potential in the QSHS is 


Aj:{z) = ByZ + Aa;(+0) = ByZ — ByX^ tanh —dN < 2 < 0 

2 Al 


( 26 ) 


where the hrst term is the vector potential of the external field, while Aa;(+0) is the solution (25) at the surface z = 0. 
Since for typical SC/QSHS structures the thickness of the QSHS ^ Al tanh (dsc/2AL), we can approximate 


A^{z) « Ha;(+0) = -ByXL tanh , -d^ < z<0. (27) 

IAl 

In the main text we often compare the energy scale associated with the magnetic field, evA^, with the superconducting 
gap A. Hence it is useful to introduce the dimensionless parameter 


^ 1^ . , „ ^_A__ 

A -B* ’ * e'uALtanh(dsc/2AL) 

where we introduced the characteristic magnetic held B* at which the 
and the flux quantum 


4’o 


(28) 


’’■^ScAl tanh (dsc/2AL) ’ 
gap closes as well as the SC coherence length 


^SC 


hv 


4>o 


h 

Ye 


(29) 


Solutions 


For solving BbcIG in the semi-inhnite plane y > 0 using the boundary conditions 


= 0 ) = ■KeiTyUyi\){y = 0 ). 


(30) 


we use the ansatz 




$fcexp 



Plugging it into the stationary Schrodinger equation we arrive at 


vhkTToT^ax + veA^-n^T^ax + - AtTr^roCTo + FTr^TxCro - A-KyTocry 




At y = 0, where (30) holds, we use 


T^zTxf^o^k = T^zTxCro{TroryCry)^k = Trz{iTz,)(Jy^k 


(31) 


(32) 


(33) 


to cast the sum {ihv/SYzTzCfy + r7r^rxCro into (F + fiv/5)T:zTx'^o- Equating this term to zero, we hnd the decay length 


hv hv 

‘ = -T = Wl- 


(34) 


The edge solution exists only for the inverted band structure with F < 0. From the boundary conditions at y = 0 we 
additionally hnd that the spinor simplihes to 


^/c = iXk, ^cTyXk, Y, io'y^fc)^ • (35) 

Now we can reformulate the problem as a system of coupled linear equations 

0 = [v{hk + eAx)ax - (y + B)(To] Xk + Aidy^k (36a) 

0 = [-v{hk + eAx)(7x - (y + E)ao] idyXk + Aiay{iay)^k (36b) 

0 = -AidyXk + [v{hk - eAx)ax + (y - B)cro] ^k (36c) 

0 = -Aiay(iay)xk + [-v{hk - eAx)(Jx + (y - E)ao] idy^k- (36d) 











Using 


i^j — “t“ 


(37) 


it is easy to show that upon multiplication by —icr^ equations (36b) and (36d) correspond to equations (36a) and 


(36c), respectively. We hence solve equations (36a I and (36c) independently. From (36a) we find 


Cfc = ^ [v{fik + eAa;)cr^ - i{pL + E)(jy\ Xk- 


(38) 


Plugging this into ( |36c[ ) we get 

0 = \{2vfik^i — 2veAxE) — e^A^) + — E^ + A^j io-y] Xk | • i^y 

0 = — {2vhkfi — 2veAxE) a^Xk + {v^{k?k^ — e^A^) + — E^ + A^) croXfc (39) 

so that Xk must be an eigenstate to Ux having the eigenvalues s = ±1 with corresponding eigenvectors (s, 1)"^. Finally 
we can solve the energy spectrum at the interface 


E^ = sveAx + a\/ + (^ — svhkY = sveAx + aWks, 

where a = ±1 and designates conduction or valence bands. Further, we defined 

Wks = -\/A2 + (/i - svhk)'^. 

For the states we obtain 

Xks — (s, 1) , '^O'yXks — (1) s) 

^ [kvka^ + (-/r =F Wks)icry] Xks 

'^^y^ks = ^ [-kvkax + (m ± hFfcs)io-o] Xks 


(40) 

(41) 


(42) 


= Jf exp 


ikx — y 


hv 


kiss’s 

'^ks4>-t 


where M comes from the normalization of the wave function. 


= (s, 1,1,-s) 


(43) 


(44) 


is the particle spinor and 


Eks 


— 1’ ^ks = ^ (^-vhk + sy± S\/A2 + (/r - svhk)'^^ 


(45) 


are the electron and hole weights, respectively. In this way we have decoupled the solution in a part coming from the 


particle solution exp 


ikx - 


0 s and a part containing the superconducting order parameter. 


The normalization consists of two parts. The first is the normalization of the particle and hole weights 

^ [WL ± (m - svHk)Wks] . 

Additionally we have to normalize over the half space y > 0 and the dimensions of the spinor 

.|rr 


= 


The full normalization is given by 


dy exp 


-2y 


hv 


, „ 9, fi'c ^hv 

• 0 = (2 + 2 s^)-—r = 2 — . 

H's \ -r y^iri |r| 


(46) 


(47) 


f^ = NN= l^Wl±{f,-shkv)WksJ^^. 


|r|' 


(48) 
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It is advantageous to rewrite the spinor (43) in the language of the coherence factors, because this form reflects the 

(49) 


natural symmetries of the system. We define the dimensionless parameter 

^ — svhk 


Vks — 


A 


which measures the deviation away from the Fermi points. With this Wks = + ? 7 ^g/|A|. The coherence factors 

are defined by (A > 0 signA = 1 |A| = A) 


1 


1 


Uks = 


|A| 


N{a,k,s) a /2 ^JWkS^ + a{fj, — svhk)Wks 
1 1 


^ \/l + vis + O^Vks^/^^ 


vis 


(50) 


Vks = 


= s 


'ks 


1 s sign A (/i — svhk + aWks) 
N{a,k,s) a/2 -^WkS^ + a{fj, — svhk)Wks 

1 Vks+aV^ + vls 

V2 

I s sign r]k, 


1 S sign 

Vks+a^l + ril^ 



iVks+a^/l+rilJ^ 


+ a\/l + C 


+ \/^ + vls-aVks^/TT^, 


In the last line we used 


1 + ??L + <^Vks\Jl + ills = (1 + vL - aVks\Jl + vls)iVks + Oi^Jl + rg^f 


(51) 


(52) 


vh 


> 0 ^ a/1 + > Iriks]■ Therefore we can simplify 


as well as that for propagating solutions (real k): riks € 1 
the sign expression to signa = a. 

Since we will be mainly interested in energy dependent quantities, like the elastic scattering and zero energy states, 
it is useful to reformulate the coherence factor in an energy dependent form. We use 


k^{E) = ^ (s/r + /3 a/(£; - seuA,)2 - A2) , /3 


= ± 1 . 


(53) 


In squaring the dispersion ( |40[ ) we loose the information a = ±, so that j3 is not in one to one correspondence with 
a. The relation can be found by checking the consistency of equation 


£1“ = sveAx + a\l A^ + (/r — svhks {E)Y = E, 
where we fix E and /?. The result is independent of /3: 

„ , , , , E > sevA^ 

® ' ' [-1, E<sevA^ 

and a is fixed by energetic constraints. In total the energy dependent formulation takes the form 


1 


^Es 


^Es 


^{E — sevA^)^ — sf5{E — sevAx)^J{E — sevA^)^ — A^ 

1 sAsign [E — sevAx] 

^ \J{E — sevAx)^ + sf3{E — sevAx)\J{E — sevA^Y — A^ 


(54) 

(55) 

(56) 

(57) 


The full wave functions are 


i’ks^y) = 


exp 


Tw 


'^ks4‘s \ ,0 / X _ 

v9 6_ E Ves^V)- 


exp 


\hV „ _ ..M 

^KesX y tiy 


N-2 


Ues^s 


(58) 
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Let us analyze the coherence factors further. For both representations it is easy to show that 


^ks \ 




Es\ 


-'Esl 


= 1 . 


( 59 ) 


Further it is interesting to note that 

= s sign [E - sevA^ , (60) 

i.e. the coherence factors are, up to a phase factor, related by /? —>■ —/3 {a —)■ —a). 

The quasi-particle character ( the particle or the hole character of the excitations) of a state can be calculated using 
the operator tt^tocto 


p0 _ 

^Es — 


(wL) - {^Es) 'I’L = sp 


Pks = (utsT uls - (VksT Ves = - 


— sevAx)'^ — 
E — sevAx 
/i — shvk 
+ (fj, — shvkY ’ 


(61) 

(62) 


showing how the quasi-particle character depends on a (/3), s and k (E). For a given energy and spin, /3 indicates the 
quasi-particle character of the state. At the Fermi points k^ = ^ {E^ = E'^p^ = sveA^ + aA) one finds 


p/? _ TDOi 

— ^k^s 


= 0 , 


(63) 


i.e. a perfect mixture of electrons and holes at the band edges. Far away from the quasi-particle gap we get 
{k = K\k\, E = e|if|, where k = e = ±1) 


,1™ Pes = ^sP, lim P^^ = Ksa, (64) 

|E|—>-oo |fc|—>-oo 


where P —>■ 1 (P —>■ — 1) indicates electrons (holes). 


DOS and sDOS 


The spin-dependent density of states in the energy interval [E, E + dE] at position y is dehned by 

Psiy,E) = -^Ti [lmG^{y,y)]^^, 

where the retarded Green’s function is given in Lehmann representation 

pfi. ,/i = V '^ksiy)i'Pksiy'))'^ 
s(y,y) 2^ E-EGk) + io+ 

k,a—:t 

- Z. (£;_i^pfc))2 + (0+)2 [E-E^{k)-^0 J. 


fe,Q; = ± 


(65) 

( 66 ) 
(67) 


The 11 component [imGf (j/, y')] corresponds to the electron block of the Green’s function. The diagonal parts of 

i’ksiy) iP’ksiy'))^ 9-^® purely real, because the y-dependence is an exponential decay instead of a plain wave. Using 
the Lorentz representation of the Dirac delta-distribution 


6{x) = lim 6eix) = lim- 

£^o+ e^0+ IT x^ + 

we find that the diagonal parts of the imaginary part of the Green’s function can be written as 


( 68 ) 


lm(Gf(y,y'))^,„ = - 0 + 

A;,a=± 


(i’ksiy) ® ii’ksiy'))^ 


(P-P«(fc))2 + (0+)2 
= E fe(2/)®(Cs(2/'))n 5{E-El{k)). 

\ / m,m 


(69) 


fc,a=± 
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Evaluating the trace we arrive at the simple form 


Ps{y.E)= Y. \m,{E){y)]^.\\^5{E-E:,{k)) 

A:,a—zb 

= E E ms{E){y)],-f5{k 

k,a=± {fco|Sj(fco)=-B} 


fco) 


<^Ei{k') 

dk' 


-1 

k'—kQ 


(70) 


Here 


bPksiE){y)]e- = ;^exp 


ikx - j/E 
hv 


f Xks 
yCTyXkt 


(71) 


is the electron part of the spinor (43). In effect it is the summation over the propagating modes at energy E 
corresponding to the quantum number s weighted with the slope of the band. The total DOS at position y is then 
given by the sum of p^{y,E) and pi{y,E). The full signal can further be obtained by integration over the half space 
y >0. 

In (70), a is chosen by the energetic constraint Eq. (55). Later on we will find that Psiy-, E) is independent of a. 

The main contribution to the spin-dependent density comes from the derivative of the band structure 


so that 


It shows poles at the band edges 


dk 


dEi[k') 


dk' 


= asvh 


-1 


svhk — p 


a/A^ -I- (/r — svhky 
1 \E — sevAx\ 


vh^(E-sevA^)^-A^' 


E = sevAx T A. 


(72) 


(73) 


(74) 


and goes to (hv) ^ away from the superconducting gap, when E — sevA^ S> A. The density is renormalized by taking 
only the j/-dependent electron contribution {Des = E — sevA^) 


r,iiE)s(E)iy) 



2 1 

o |rr 

e~ 


^ hv_ 


= exp 

= exp 




(s^-b 1-f 1-b s^) 


Ea 2 

hv 


1 


Ea 2 

hv 


Dl^ + aP^Dl^-A^\DEs 

1 

Dl^+dy^DY^DEs' 


(75) 


Again the a dependence is obsolete, since - accordingly to (55) - \Des\ = \E — sevA^l cancels the choice of a. Hence 

-1 


ps{y.E) = Y iKf(£)(y) 


3=± 


dEf{k') 


dk' 


k'=kP{E) 


(0[-D - A]-b 0[D - A]). 


(76) 


The Heavyside-0 functions ensure that E does not he within the spin-dependent gap (74), which was previously 
enforced by the Dirac-(5 distribution. 

Since (73) is independent of /3, the summation simplifies to a factor of two. In total we find the closed form 

\E — sevAA 


2 |r| 

Ps{y,E) =-^exp 




1 


tiv \ hv ^[E - sevAa^Y - A^ 
Integration over the half space gives 

\E — sevA^-l 


{Q[sevAx — A — E] -b 0[if — {sevA^ -b A)]). (77) 




— sevA^)"^ — A^ 


{Q[sevAx — A — E] -b 0[if — (sevA^ + A)]). 


(78) 
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It is easy to see that the s-polarization Pg = — pi vanishes for = 0, because all s dependence vanishes. But also 

for A = 0 we find vanishing since 


1 \E — sevA^\ 
fiv yj{E — sevAx)'^ — 


1 

— = const. 
nv 


(79) 


together with Q[E — sevA^] + Q[sevA^ — E] = 1, which is expected for an ungaped system. Again p^ and pi are 
equivalent and Ps is zero. 

Interestingly the spin-dependent density is independent of the chemical potential p. Comparing with the dispersion 
(401 we find that p just spin-dependently shifts k, while the density of states is calculated by summing over all k values 
at a certain energy. At zero vector potential the system shows to be a gaped superconductor, cf. Fig. Increasing 


Ax closes this gap. In Figj^we analyze the gap closing in more detail. For zero vector potential (cf. Figs. 6a and 6b) 
we find that states of a = -I- and a = — { conduction and valence bands) are separated by a superconducting gap. 
Additionally we find a Dirac cone between states of same a but opposite s. The cone is protected because states of 
different s are orthogonal. Increasing A^ shifts the branches of different s in opposite directions and therefore the gap 


opens around energies E = sevAx, while k = Epj(hv) is fixed. This is indicated by the vertical grid lines, cf. Fig. 6c 


Increasing the vector potential (see Fig. 6d) is reflected by shifting part of the density into the superconducting gap. 
This is accompanied by a spin-dependent signal, since the states of different s do no longer live at the same energy. 
When |enAa;| = |A|, i.e. \By\ = in Figs. 6e andIM the indirect gap in the dispersion is closed, indicated by a peak 


in the DOS at zero energy. However, the spin branches have a gap. Increasing the vector potential even further shifts 
the relative spin branches far enough that they switch energetic order at By 
and Uhl 


2.255* for our parameters, cf. Figs 6g 


Transport within the BTK formalism 


A way to measure the spin polarization of the edge states is transport through an interface between a QSHS and 
the QSHS in the proximity to the s-wave superconductor (NSC-junction). A possible realization of such a junction is 
the Y-forked NSC junction described in the main text. Due to spin-momentum locking the incoming electronic state 
on the left side of the sample has a fixed spin on a certain edge of the sample. The normal part of the junction is 

described by Pl^dG with A —>• 0. Using \/D^^ + IS? ^ \Des\ (^1 — 2 % ) 


lim = 
A-s-O 


A 


Es 


/3 ssignDes 

lim = - 

A-s-O 


_- < ’ sign[i:)£)s]s/3 < 0 

\JD%^{1 - s,dsign[D£;s]) + s/SsigniDBs]^ sign[i:)£)s]s/3 > 0 

A 


+ s/?sign[DBs]) - sl3sign[DEs 


I A2 


_ I 0 , sign[D£;s]s ,5 > 0 

[ssign[i:)£;s], sign s/3 < 0 
The solutions in the normal regime hence take the form 

electrons: '0^/ \y) 

holes: 


exp 


^ y hv 


exp 


iKes X y 


ssign[Ds„](/)_, 


where /3 no longer is a degree of freedom but fixed by (80). For incoming electrons 


/3(e ) = s sign[DEs] 


(80) 


(81) 


(82) 


has to be chosen. The helical nature of the edge states forbids electron reflection as long as the tC symmetry is 
preserved. The only possible scattering channels are the electron transmission and reflection as a hole (Andreev 
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jU=2A, By=OB^ 




E[A] 


(a) E(k), By = OS* 


(b) p(E), By = OB, 


jU=2A, By=0.5B^ 




k [^sc] E [A] 


(c) E(k), By = 0.5S* 


(d) p(E), By = 0.5B* 


jU=2A, By=lB^ 




E[A] 


(e) E(k), By = IS* 


(f) piE), By = IS* 


fi=2A, By=2.l5B^ 



jU= 2 A, £3 =2.155* 



E[A] 


(g) E{k), By = 2.15S* 


(h) p(E), By = 2.15S* 


FIG. 6: The gap closing is analyzed by comparing the dispersion (left column) with the DOS (right column) for different values 
of By. In the left column s = 1 (—1) is indicated by solid (dashed) curves. The particle character is color encoded, blue (red) 
being used for hole (electron) like states. 
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reflection). When the excitation energy lies within the superconducting gap electron transmission must be zero and 
perfect Andreev reflection was predicted for both geometries, provided that the width of the ribbon is large enough to 
separate the counter propagating edge channels [23 EH] ■ From matching the wave functions at the interface at a; = 0 


+ m 


0 

s sigTi[DEs](l)-, 


= T 


3 

Es 


^Es^- 


(83) 


we find 


RUE) = s sign)!!^;,,]^. 


(84) 


^Es 


where R\{E) and T are the Andreev reflection and electron transmission amplitudes, respectively. We investigate 
the spin-dependent currents within the BTK formalism |2S|. When the voltage V is applied across the junction, the 
current can be calculated by 


OO 

^btk = ^ J dA(l + |i?i|2)(/(A-eF)-/(A)). 


(85) 


Here f{E) is the distribution function of the carriers at excitation energy E and temperature T 

1 


The excess current is defined by 


f{E) = 


exp 


E 

ksT 


1 


( 86 ) 


E = - 
h 


OO 

J dE\R^^mf{E-eV)-fiE)) 


(87) 


It gives the current which is transported additionally to the ohmic one and is generally measured at very high voltages. 
In the main text we found that the spin-dependent excess current shows a maximum at hnite By. To analyze this 
further let us have a look at the excess current at zero temperature 


eV 


eV+sgtlisBy 


T=0 e 


= - dE\R%\^ = - 


h 


h 


dDEs\RA\ , 


( 88 ) 


+sg,iJ.^By 


where in the last step we used that in the Andreev reflection coefficient E and By only appear via Des = E + sg:tfJ-BBy. 
g* is the effective g-factor defined in the main text and the Bohr magneton. With 

|^^|2 ^ Des- 3j3{e-)^/Dl^- ^ Des - signJi;,^]^ - ^ 


(90) 


DEs-syie-)^D'i^-A^ 


Des - sigaDEs\/- A^ 

DEs + sP{e-)^Dl^-A^ 


Des + signDEs\/Dy^ - 


the maximum in By can be found by setting 

dsyle^ oc \R\{eV + S 5 */rBHy)p - \R\{+sg^g,BBy)\'^ = 0. 
With \RW’^{Des) = \Ra\’^{—Ees) we find that this equation can be solved by 

nmax _ „ 

" 2g,^XB ■ 


Another typical way to measure non-local conductances are d//dF characteristics. From (85) we obtain 

OO 

(df{E-eVy 


dr 

dV 


dE{l + \Ry\y 




dV 


(91) 


(92) 




















